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Upper limits on the spin-independent (SI) as well as spin-dependent (SD) elastic scattering cross 
sections of low mass (~ 2 - 20 GeV) WIMPs (Weakly Interacting Massive Particles) with protons, 
imposed by the upper limit on the neutrino flux from WIMP annihilation in the Sun given by the 
Super-Kamiokande (S-K) experiment, and their compatibility with the "DAMA-compatible" regions 
of the WIMP parameter space — the regions of the WIMP mass versus cross section parameter space 
within which the annual modulation signal observed by the DAMA/LIBRA experiment is compatible 
with the null results of other direct detection experiments — are studied within the frame work of 
a self-consistent model of the finite-size dark matter (DM) halo of the Galaxy. The halo model 
includes the gravitational influence of the observed visible matter of the Galaxy on the phase space 
distribution function of the WIMPs constituting the Galaxy's DM halo in a self-consistent manner. 
Unlike in the "Standard Halo Model" (SHM) used in earlier analyses, the velocity distribution of the 
WIMPs in our model is non-Maxwellian, with a high-velocity cutoff determined self-consistently by 
the model itself. The parameters of the model are determined from a fit to the rotation curve data of 
the Galaxy. We find that, for our best fit halo model, for SI interaction, while the S-K upper limits 
do not place additional restrictions on the DAMA-compatible region of the WIMP parameter space 
if the WIMPs annihilate dominantly to bb and/or cc, portions of the DAMA-compatible region 
can be excluded if WIMP annihilations to r"'"r~ and vu occur at larger than 35% and 0.4% levels, 
respectively. For SD interaction, on the other hand, the restrictions on the possible annihilation 
(— I I channels are much more stringent: they rule out the entire DAM A region if WIMPs annihilate to 

^ J t'^t~ and ui? final states at greater than ~ 0.05% and 0.0005% levels, respectively, and/or to 66 and 

cc at greater than ^-^0.5% levels. The very latest results from the S-K Collaboration [T. Tanaka et 
al, Astrophys. J. 742:78 (2011)] make the above constraints on the branching fractions of various 
WIMP annihilation channels even more stringent by roughly a factor of 3-4. 
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I. INTRODUCTION 

Weakly Interacting Massive Particles (WIMPs) (hereafter generically denoted by x) with masses m^^ in the range of 
few GeV to few TeV are a natural candidate for the dark matter (DM) in the Universe; See e.g., Refs. [1-5] for reviews. 
Several experiments are currently engaged in efforts to directly detect such WIMPs by observing nuclear recoils due to 
scattering of WIMPs off nuclei in suitably chosen detector materials in underground laboratories. Recent results from 
some of these direct detection (DD) experiments, in particular the annual modulation of the nuclear recoil event rates 
reported by the DAM A/LIBRA collaboration [6] and the excess of low energy recoil events reported by the CoGeNT 
collaboration [7] have raised the interesting possibility [8, 9] that these events could be due to WIMPs of relatively 
low mass, approximately in the range ~ 5-10 GeV, interacting with nuclei with a WIMP-nucleon spin- independent 
elastic cross section in the region of few x lO"** pb, without conflicting with the null results from other experiments 
such as XENONIO [10], XENONIOO [11] and CDMS-II-Si [12]. Earlier analyses (before the CoGeNT results [7]) had 
also found similar compatibility of the DAMA/LIBRA annual modulation signal with the null results from other DD 
experiments; sec, e.g., Refs. [13-15]. ^ 

Scattering of WIMPs off nuclei can also lead to capture of the WIMPs by massive astrophysical bodies such as the Sun 
or the Earth if, after scattering off a nucleus inside the body, the velocity of the WIMP falls below the escape velocity 
of the body. The WIMPs so captured over the lifetime of the capturing body would gradually settle down to the core 
of the body where they would annihilate and produce standard model particles, e.g., W~^W~ , Z^Z^ ^ t~^t~ , ti, bb, cc, 
etc. Decays of these particles would then produce neutrinos, gamma rays, electrons-positrons, protons-antiprotons, 
etc. For astrophysical objects like the Sun or the Earth, only the neutrinos would be able to escape. Detection of these 
neutrinos by large neutrino detectors can, albeit indirectly, provide a signature of WIMPs. Although no detection 
has yet been reported, the Super-Kamiokande (S-K) detector, for example, has provided upper limits on the possible 
neutrino flux from WIMP annihilation in the Sun as a function of the WIMP mass [18-20]. Similarly, the 7-rays 
produced in the annihilation of the WIMPs in suitable astrophysical environments with enhanced DM density but 
low optical depth to gamma rays, such as in the central region of our Galaxy, in dark matter dominated objects such 
as dwarf galaxies, and in clusters of galaxies, can offer a complimentary avenue of indirect detection (ID) of WIMPs. 
Although no unambiguous gamma ray signals of dark matter origin have been reported, a recent analysis [21] of 
the spectral and morphological features of the gamma ray emission from the inner Galactic Center region (within a 
Galactocentric radius of ~ 175 pc) measured by the Fermi Gamma-ray Space Telescope (FGST) seems to suggest the 
presence of a gamma ray emission component which is difficult to explain in terms of known sources and/or process 
of gamma ray production, but is consistent with that expected from annihilations of WIMPs of mass in the 7-9 GeV 
range (annihilating primarily to tau leptons) with a suitably chosen density and distribution of the dark matter in 
the Galactic Center region; see, however, Ref. [22] for a different view. 

In this paper we focus on the neutrinos produced by annihilations of WIMPs in the core of the Sun, and study the 
constraints imposed on the WIMP mass vs. WIMP-nuclcon cross section, for low-mass ( 20 GeV) WIMPs, from 
non-detection of such neutrinos. This is done within the context of a self-consistent model of the finite-size dark halo of 
the Galaxy [15, 23] that includes the gravitational effect of the observed visible matter on the DM in a self-consistent 
manner, with the parameters of the model determined from fits to the rotation curve data of the Galaxy [24, 25]. 

The expected flux of neutrinos from the Sun due to WIMP annihilations depends on the rate at which WIMPs are 
captured by the Sun. The capture rate depends on the density as well as the velocity distribution of the WIMPs 
in the solar neighborhood as the Sun goes around the Galaxy. The density and velocity distribution of the WIMPs 
in the Galaxy are a priori unknown. Most earlier studies of neutrinos from WIMP capture and annihilation in the 
Sun have been done within the context of the so-called "Standard Halo Model" (SHM) in which the DM halo of the 
Galaxy is described by a single component isothermal sphere [26] with a Maxwellian velocity distribution of the DM 
particles in the Galactic rest frame [1, 27, 28]). The velocity distribution is isotropic, and is usually truncated at a 
chosen value of the escape speed of the Galaxy. The density of DM in the solar neighborhood is typically taken to be 
in the range Pdm,© 0.3 ± 0.1 GeV/ cm^ [29-32] ^. The velocity dispersion, (u^)^/^. the parameter characterizing the 
Maxwellian velocity distribution of the SHM, is typically taken to be ~ 270kms~^. This follows from the relation [26], 



^ The question of compatibility of the DAMA/LIBRA and CoGeNT results with the null results of other experiments, however, remains 
controversial; see, e.g., the results of a recent reanalysis of the CDMS-II Germanium data with a lowered recoil-energy threshold of 2 
keV [16], as well as the recent results from the XENONIOO collaboration [17], both of which claim to disfavor such a compatibility. 

^ See, however, recent analyses [33, 34] which claim a value closer to 0.4 GeV/ cm'^ with uncertainty ^ 10%. 
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(„2^1/2 ^ ^3^^^^^ 

between the velocity dispersion of the particles constituting a single-component self-gravitating 

isothermal sphere and the asymptotic value of the circular rotation speed, Uc,oo, of a test particle in the gravitational 
field of the isothermal sphere and assuming Wc,oo ~ Wc,© ~ 220 km s""'^, where Vc,q is the measured value of the circular 
rotation velocity of the Galaxy in the solar neighborhood. Neutrino flux from DM annihilation in the Sun for 
low mass WIMPs and the resulting constraints on WIMP properties from the Super-Kamiokande upper limits on 
such neutrinos have been studied within the context of the SHM in Refs. [19, 36-38], which showed that the Super- 
Kamiokande upper limits on the possible flux of neutrinos from the Sun place stringent restrictions on the DAMA 
region of the WIMP parameter space. 

Whereas the SHM serves as a useful benchmark model, there are a number of reasons why the SHM does not provide 
a satisfactory description of the dynamics of the Galaxy. First, it does not take into account the modification of the 
phase space structure of the DM halo due to the significant gravitational effect of the observed visible matter on the 
DM particles inside and up to the solar circle. Second, the isothermal sphere model of the halo is infinite in extent and 
has a formally divergent mass, with mass inside a radius r, M{r) oc r, as r — > oo, and is thus unsuitable for representing 
a halo of finite size. Third, the procedure of truncating the Maxwellian speed distribution at a chosen value of the 
local (solar neighborhood) escape speed is not a self-consistent one because the resulting speed distribution is not 
in general a self-consistent solution of the steady-state collisionless Boltzmann equation describing a finite system 
of collisionless DM particles. In addition, since the rotation curve for such a truncated Maxwellian is, in general, 

not asymptotically fiat, the relation (w^)^/-^ = \J\vc^oo used to determine the value of (w^)^/^ in the Maxwellian 
speed distribution of the isothermal sphere, as done in the SHM, is not valid in general. Finally, recent numerical 
simulations [32] seem to find that the velocity distribution of the Dark Matter particles deviates significantly from the 
usual Maxwellian form. These issues are further discussed in detail in Rcf. [15], where we discussed a self-consistent 
model of the finite-size dark halo of the Galaxy that avoids the above mentioned inconsistencies of the SHM and also 
studied the constraints on WIMP properties from the results of the direct detection (DD) experiments within the 
context of this self-consistent halo model. It is of interest to extend this study to the case of indirect detection (ID) 
of WIMPs via neutrinos from WIMP annihilations in the Sun, which is the purpose of this paper. 

Our model of the phase space structure of the finite-size DM halo of the Galaxy is based on the so-called "lowered" 
(or truncated) isothermal models (often called "King models") [26] of the phase-space distribution function (DF) 
of collisionless particles. These models are proper self-consistent solutions of the collisionless Boltzmann equation 
representing nearly isothermal systems of finite physical size and mass. There are two important features of these 
models: First, at every location within the system a DM particle can have speeds up to a maximum speed which is 
self-consistently determined by the model itself. A particle of maximum velocity at any location within the system 
can just reach its outer boundary, fixed by the truncation radius, a parameter of the model, where the DM density 
by construction vanishes. Second, the speed distribution of the particles constituting the system is non-Maxwellian. 
To include the gravitational effect of the observed visible matter on the DM particles, we modify the "pure" King 
model DF by replacing the gravitational potential appearing in the King model DF by the total gravitational potential 
consisting of the sum of those due to DM and the observed visible matter. This interaction with the visible matter 
influences both the density profile and the velocity distribution of the dark matter particles as compared to those for 
a "pure" King model. In particular, the dark matter is pulled in by the visible matter, thereby increasing its central 
density significantly. When the visible matter density is set to zero and the truncation radius is set to infinity, our 
halo model becomes identical to that of a single-component isothermal sphere used in the SHM. For further discussion 
of the model, see [15, 23]. 

The DM distribution in the Galaxy may have significant amount of substructures which may have interesting effects 
on the WIMP capture and annihilation rates [39]. However, not much information, based on observational data, is 
available about the spatial distribution and internal structures of these substructures. As such, in this paper we shall 
be concerned only with the smooth component of the DM distribution in the Galaxy described by our self-consistent 
model mentioned above, the parameters of which are determined from the observed rotation curve data for the Galaxy. 

The non-Maxwellian nature of the WIMP speed distribution in our halo model makes the calculation of the WIMP 



^ A somewhat higher value of fc,0 ~ 250 km s ^, as suggested by a recent study [35], would imply a correspondingly higher value of 
(t)2)i/2._^^ ~ 306kms-i. 
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capture (and consequently annihilation) rate non-trivial since the standard analytical formula for the capture rate given 
by Gould [40] and Press and Spergel [41], which is widely used in the literature, is not valid for the non-Maxwellian 
speed distribution in our halo model, and as such has to be calculated ab initio; see section III. 

We calculate the 90% C.L. upper limits on the WIMP-proton spin-independent (SI) as well as spin-dependent (SD) 
elastic cross sections as a function of the WIMP mass, for various WIMP annihilation channels, using the 90% C.L. 
upper limits on the rates of upward-going muon events due to neutrinos from Sun derived from the results of S-K 
collaboration ([18], [19] and references therein). * We then study the consistency of those limits with the 90% C.L. 
"DAMA-compatible" regions — the regions of the WIMP mass versus cross section parameter space within which 
the annual modulation signal observed by the DAM A/LIBRA experiment [6] is compatible with the null results of 
other DD experiments — determined within the context of our halo model [15]. We find that the requirement of such 
consistency imposes stringent restrictions on the branching fractions of the various WIMP annihilation channels. For 
example, in the case of spin-independent WIMP-proton interaction, while the S-K upper limits do not place additional 
restrictions on the DAMA-compatible region of the WIMP parameter space if the WIMPs annihilate dominantly to 
bb, cc, portions of the DAMA-compatible region can be excluded if WIMP annihilations to r+r^ and fD occur at 
larger than 35% and 0.4% levels, respectively. In the case of spin-dependent interactions, on the other hand, the 
restrictions on the branching fractions of various annihilation channels are much more stringent. Specifically, they 
rule out the entire DAMA region if WIMPs annihilate to r+r" and lyp final states at greater than ~ 0.05% and 
0.0005% levels, respectively, and/or to bb and cc at greater than ^0.5% levels. ^ The very latest results from the 
S-K Collaboration [20] make the above constraints on the branching fractions of various WIMP annihilation channels 
even more stringent by roughly a factor of 3-4. 

The rest of the paper is organized as follows: In section II we briefly describe the self-consistent model of the DM halo 
of the Galaxy. The formalism of calculating the WIMP capture and annihilation rates in the Sun within the context 
of our DM halo model, and that for calculating the resulting neutrino flux and event rate in the Supcr-Kamiokande 
detector, are discussed in sections III and IV, respectively. Our results and the constraints on the WIMP properties 
implied by these results are described in section V. The paper ends with a Summary in section VI. 

II. THE SELF-CONSISTENT TRUNCATED ISOTHERMAL MODEL OF THE MILKY WAY'S DARK 

MATTER HALO 

The phase space distribution function (DF) of the DM particles constituting a truncated isothermal halo of the Galaxy 
can be taken, in the rest frame of the Galaxy, to be of the "King model" form [15, 23, 26], 

/(x,v) . m = I p^(^--'r'^' {-'^"^ - 1) ^ > 0' (1) 

[ for f < , 

where 

£:(x) ^ ^n) + $(x)) ^ vl;(x) - , (2) 

is the so-called "relative energy" and ^(x) = — $(x) -|-$(rt) the "relative potential" , $(x) being the total gravitational 
potential under which the particles move, with boundary condition $(0) = 0. The relative potential and relative 
energy, by construction, vanish at ]x] = rt, the truncation radius, which represents the outer edge of the system where 
the particle density vanishes. At any location x the maximum speed a particle of the system can have is 

l'max(x) = V2«'(x) , (3) 



* After the completion of the main calculations of the present work, new results of the S-K collaboration's search for upward-going muons 
due to neutrinos from Sun [20] have appeared. We include, at the end of section V, a discussion of the new results of Ref. [20] and the 
resulting constraints on various WIMP annihilation channels. 

^ In the present paper, the CoGeNT results [7] are not included in the analysis. Preliminary results of the analysis [42] to find the 
"CoGeNT-compatible" region in the WIMP mass vs. cross section plane within the context of our halo model indicates that its inclusion 
will not significantly change the above constraints on the branching fractions for the various annihilation channels. 
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at which the relative energy 8 and, as a consequence, the DF (1), vanish. The model has three parameters, namely, 
Pi , a and . Note that the parameter tr in the King model is not same as the usual velocity dispersion parameter of 
the isothermal phase space DF [26]. Also, in our calculations below, we shall use the parameter Pdm,0, the value of 
the DM density at the location of the Sun, in place of the parameter pi . 



Integration of /(x, v) over all velocities gives the DM density at the position x: 



Pom(x) = 



Pi 



2^3/2 



Pi 



which satisfies the Poisson equation 



(27rCT2) 

exp 



v'2*(x) 



VI/ (x) 



erf 



exp 



*(x) - 



4*(x) 



- 1 

2*(x 
3o-2 



V2$om(x) =47rG/9oM(x), 



(4) 
(5) 

(6) 



where $dm is the contribution of the DM component to the total gravitational potential, 

$(x) = $om(x) + $vm(x), (7) 

in presence of the visible matter (VM) whose gravitational potential, $vm, satisfies its own Poisson equation, namely, 

V2$vm(x) = 47rGpvM(x) . (8) 

We choose the boundary conditions 



$dm(0)-$vm(0) = 0, and (V$om)|x|^o = (V* 



vm;|x|=o 



0. 



(9) 



The mass of the system, defined as the total mass contained within rt, is given by GM{rt)/rt — [$(00) — ^{rt)]. 
Note that, because of the chosen boundary condition $(0) = 0, $(00) is a non-zero positive constant. 

Since the visible matter distribution Pvm(x), and hence the potential $vm(x), are known from observations and 
modeling, the solutions of equation (6) together with equations (5), (7) and the boundary conditions (9), give us a 
three-parameter family of self-consistent pairs of Pom(x) and $om(x) for chosen values of the parameters (pi ,cr,rt). 
The values of these parameters for the Galaxy can be determined by comparing the theoretically calculated rotation 
curve, Vc{R), given by 



vl{R)^R—[<fiR,z = 0) 



R 



d_ 

dR 



$oM(i?,^-0) + $vM(i?,2-0) 



(10) 



with the observed rotation curve data of the Galaxy. (Here R is Galactocentric distance on the equatorial plane and 
z is the distance normal to the equatorial plane.) This procedure was described in detail in Refs. [15, 23] where, 
for the visible matter density distribution described there, we determined the values of the parameters and a 
that gave reasonably good fit to the rotation curve data of the Galaxy [24, 25] for each of the three chosen values 
of the parameter Pdm,© = 0.2, 0.3 and 0.4 GeV/cm'^. These models are summarized in Table I, which we use for 
our calculations in this paper. The density profiles, mass profiles, velocity distributions of the DM particles and the 



Model 


Pdm,o 


rt 


a 




(GeV/cm^) 


(kpc) 


(kms^^) 


Ml 


0.2 


120.0 


300.0 


M2 


0.3 


80.0 


400.0 


M3 


0.4 


80.0 


300.0 



TABLE I: Parameters of our self-consistent model of the Milky Way's Dark Matter halo that give good fits to the Galaxy's 
rotation curve data, for the three chosen values of the DM density at the solar neighborhood 



resulting rotation curves in each of these models are discussed in detail in Ref. [15]. 



With our halo model specified, we now briefly review the basic formalism of calculating the WIMP capture and 
annihilation rates within the context of our halo model. 
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III. CAPTURE AND ANNIHILATION RATES 



The capture rate per unit volume at radius r inside the Sun can be written as [40, 41] 

^(r)= Jd^uMn^n-iw), (11) 

where /(u) is the WIMP velocity distribution, as measured in the Sun's rest frame, in the neighborhood of the 
Sun's location in the Galaxy, and w{r) = ^/vF^\^w^^Jr) is the WIMP's speed at the radius r inside the Sun, 
'Wosc(?') being the escape speed at that radius inside the Sun, which is related to the escape speed at the Sun's core, 
Wosc.corc ~ 1354 kms~^, and that at its surface, Wesc.smf ~ 795kms~^, by the approximate relation 

M(r) 

'^csc (^) — (^csc,corc) [(^csc,coro) ('^csc,surf) ] ■ (-^^) 

The quantity n~{w) is the capture probability per unit time, which is just the product of the scattering rate and the 
conditional probability that after a scattering the WIMP's speed falls below the escape speed. 



We shall here consider only the elastic scattering of the WIMPs off nuclei. The dominant contribution to the WIMP 
capture rate will come from the WIMPs scattering off hydrogen and helium nuclei. While for hydrogen, both spin- 
independent (SI) as well as spin-dependent (SD) cross sections, cr^p and cr^p , respectively, will contribute, only SI 
cross section for helium is relevant. (We neglect here the small contribution from '^He). In general, the effective 
momentum-transfer (q) dependent WIMP-nuclcus SI scattering cross section, <T^^{q), can be written in the usual way 
in terms of the "zero- momentum" WIMP-proton (or WIMP-neutron) effective cross section, cr^p = cr^Jj, as 



id) 



(13) 



where A is the number of neutrons plus protons in the nucleus, /x^^a and /i^p are the reduced masses of WIMP-nucleus 



and WIMP-proton systems, respectively, with = {rrii my^)/{mi 
Fifi) = 0) which can be chosen to be of the form [1] 



), and F{q ) is the nuclear form- factor (with 



exp 



exp 



(14) 



Here R 



0.91 



rUA 



1/3 



0.3 



GeV/ 

nuclear coherence energy, tua being the mass of the nucleus 



X 10 ^ cm is the nuclear radius and Eq = 3h /{2mAR ) is the characteristic 



With the above form of the nuclear form factor, the kinematics of the capture process [40] allows us to write the 
capture probability per unit time, n~{w), as 



n- 



UA O 



XA 2£'o 



w 



exp 



2En 



exp 



2E, 







,,2 



(15) 
/i± 1 



jn 

where ua is the number density of the scattering nuclei at the radius r inside the Sun, and /i = — — , ^± 

niA ^ 

The 9 function ensures that those particles which do not lose sufficient amount of energy to be captured are excluded. 



We shall use Equation (15) to calculate {w) for helium (A = 4). For hydrogen, however, there is no form-factor 
suppression, and the expression for fl~('w) is simpler: 

Hydrogen : (w) = (w^^ - ^uA 9 (wl^ - ^uA , (16) 

where tih is the density of hydrogen (proton) at the radius r inside the Sun. Note that in equations (15) and (16), 
the quantities w, Wcsc, n-A and hh are functions of r. 



7 



The WIMP velocity distribution appearing in equation (11) is related to the phase space DF defined in equation (1) 
(valid in the rest frame of Galaxy) by the Galilean transformation 

/(u) = — /(x = X0,v = u + V0) , (17) 

where x© represents the sun's position in the Galaxy [R = 8.5kpc,z = 0) and is the Sun's velocity vector in the 
Galaxy's rest frame. Note that Gould's original calculations and the final formula for the WIMP capture rate given in 
Ref. [40] , which are widely used in the literature, use a Maxwellian velocity distribution of the WIMPs in the Galaxy, 
and as such, cannot directly be used here since the WIMP velocity distribution in our case is non-Maxwellian. In 
particular, note that the DF / of equation (1) vanishes for speeds v > Wmax defined in equation (3). Consequently, 
equation (11) above can be written as 

dC 27r ;-timax(cos6() 

— — (r) = / d(cos0) / u dw/ (x = xq, V = u + vq) wfi^(u;) , (18) 

dV m^J-i A„„„=^0 

where vq « 220 - 250kms~^ is the Sun's circular speed in the Galaxy, and Umax is given by the positive root of the 
quadratic equation 

"max + 2Umax Vq COS 6* = 2^'(x = Xq) . (19) 

The total WIMP capture rate by the Sun, C©, is given by 

Co= ["%nr'dr^, (20) 







dV 



where Rq is the radius of the Sun. 



In this work we shall neglect the effect of evaporation of the captured WIMPs from the Sun ^, and make the standard 
assumption that the capture and annihilation processes have reached an approximate equilibrium state over the long 
lifetime of the solar system (t© ^ 4.2 billion yrs). Under these assumptions, the total annihilation rate of WIMPs in 
the Sun is simply related to the total capture rate by the relation 

To « (21) 

IV. NEUTRINO FLUX FROM WIMP ANNIHILATION IN THE SUN AND EVENT RATE IN THE 

DETECTOR 

A. The neutrino energy spectrum 

In this subsection we collect together the known results for the energy spectra of neutrinos emerging from the Sun, 
for various WIMP annihilation channels [43, 44, 46, 47], for use in the calculations described in the next subsection. 

The differential flux of muon neutrinos observed at Earth is [43] 

d^^\ ^ fdNA _ 

rf^j-i^E^4^j^' = (22) 

where Tq is the rate of WIMP annihilation in the Sun, D is the Earth-Sun distance, F stands for the possible 
annihilation channels. Bp is the branching ratio for the annihilation channel F and (^j^^ is the differential 
energy spectrum of the neutrinos of type i emerging from the Sun resulting from the particles of annihilation 



Note, however, that evaporation may not be neghgible for WIMP masses below ~ 4 GeV depending on the magnitude of the annihilation 
cross section [36]. 
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channel F injected at the core of the Sun. WIMPs can annihilate to all possible standard model particles e.g. 
e+e", /i^/i^, T+T^, Vei'e, t^/j.i'fj,, ^t-Vt, qq pairs and also gauge and higgs boson pairs (M^+M^", ZZ, hh), etc. In this 
paper we are only interested in low mass (^2-20 GeV) WIMPs. Therefore, we will not consider WIMP annihilations 
to higgs and gauge boson pairs and top quark pairs. Light quarks like u, d, s contribute very little to the energetic 
neutrino flux [44], and are not considered. The same is true for muons. So, in this paper we consider only the channels 
T+r" , bb, cc and Pv. 

The neutrino energy spectra, ^^^^ , have been calculated numerically (see, e.g., [44, 45]) by considering all the 

details of hadronization of quarks, energy loss of the resulting heavy hadrons, neutrino oscillation effects, neutrino 
energy loss due to neutral current interactions and absorption due to charged current interactions with the solar 
medium, i/r regeneration, etc. However, the numerical results in [44, 45] are given for WIMP masses to^ > 10 GeV, 
and it is not obvious if those are valid for lower WIMP masses which are of our primary interest in this paper. In 
any case, given the presence of other uncertainties in the problem, particularly those associated with astrophysical 
quantities such as the local density of dark matter and its velocity distribution, we argue that it is good enough 
to use — as wc do in this paper — approximate analytical expressions for the neutrino spectra available in the 
literature [43, 46, 47]. We are interested in the fluxes of muon neutrinos and antineutrinos, for which we use the 
analytic expressions given in Ref. [43], which neglect neutrino oscillation effects. By comparing with the neutrino 
fluxes obtained from detailed numerical calculations [44] , we find that for small WIMP masses below ^ 20 GeV (the 
masses of our interest in this paper), the analytic expressions for the muon neutrino fluxes given in [43] match with 
the results of detailed numerical calculations [44] to within a few percent. 

The main effect of the interaction of the neutrinos with the solar medium is that [47] a neutrino of type i (= z/^, i/^) 
injected at the solar core with energy El°'° emerges from the Sun with an energy Ei given by 

Er ^Ej{l-E,n), (23) 

and with probability 

= (1 + Ern)-"' = (1 - Emr , (24) 

with 

a^^ 5.1 , ap^ ^ 9.0 , t^^ = 1.01 x 10"'"^ GeV"^ , and t^^ = 3.8 x 10"^ GeV"^ . (25) 

Below we write down the expressions for the energy spectra of neutrinos emerging from the Sun for the four annihilation 
channels considered in this paper: 



1. T^T channel : Neutrinos from decay of r leptons (r — i- ^v^j^v^) 

For this channel, the spectrum of muon-typc neutrinos at the solar surface, including the propagation effects in the 
solar medium, can be written as [43] 

where the relationship between Ei and El"" ^ and the values of on and r^, are as given by equations (23) and (25), 
respectively, and 

is the neutrino spectrum due to decay of the r-lcptons injected at the solar core by WIMP annihilations. Here 
Tt-^i_lu^u^ — 0.18 , and E± = — j-^ — ^ with 7 = (l — /3^) =^ m^/nir , ttLt being the r-lcpton mass. 



Note that the j^t-s produced from r decay may again produce rs by charged current interactions in the solar medium, 
and these secondary rs can decay to give secondary v^s. But these i^^s would be of much lower energy compared to 
the primary v^s from r decay, and are not considered here. 
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2. bb channel : Neutrinos from decay of &-quark hadrons (6 c^v^) 

The treatment is similar to the case of t decay described above. However, here the hadronizatioii of quarks and 
stopping of heavy hadrons in the solar medium have to be taken into account. The resulting spectrum of muon- 
neutrinos emerging from the Sun is given by [43] 

JAT \ rEa / 1 rlN \^^'^^°''^ /'rIM''°'°\ 



where m;, is the 6-quark mass, « 0.71m-^ is the initial energy of the b-quark hadron (the fragmentation function 
is assumed to be a sharply peaked function [43]), 

is the neutrino spectrum resulting from decay of the 6-quark hadron injected at the solar core, and 



hadron 7-^ r j-^ j-^ 



1 rfiV\"7'"" , E, 



Eo Ed 



(30) 



with Ec ~ 470 GeV, is the distribution of the hadron's energy at the time of its decay if it is produced with an initial 
energy Eq. In equation (29), Th^fj_^^x ~ 0.103 is the branching ratio for inclusive semi-lcptonic decay of 6-quark 

hadrons to muons [48], and E± = — ^ with 7 = (l — /3^) = Ed/irn, . 

7(l=F/3) 

3. cc channel : Neutrinos from decay of c-quark hadrons (c s^v^) 

Again, this is similar to the case of 6-decay discussed above, except that the kinematics of the process is slightly 
different. The resulting muon neutrino spectrum is given by [43] 

where TOc is the c-quark mass, Eq « 0.55m^ is the initial energy of the charmed hadron. 



is the neutrino spectrum resulting from decay of the c-quarks injected at the solar core, with Fc^^i/^x = 0.13, 

E± = , 7 = (1 - P^'^^ = Ed/m, , and - — I (Sq, Ed) is given by equation (30) with E, « 250 GeV 

7(1=FP) \NdEdJ 
for c-quark. 

4. vv channel : (xx i'm'^m) 



In this case the spectrum of muon neutrinos emerging from the Sun is simply given by 



= {1-E^ny , (i = Vt,,Vf,) (33) 



where 

'dN. 
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B. Calculation of Event Rates in the Super-Kamiokande detector 

The rate of neutrino induced upward-going muon events, TZ, in the S-K detector due to i^^s and P^s from WIMP 
annihilation in the Sun can be written as 

E 1 1 ^^^iE.,y)VME,)rf;^"dE,dy, (35) 

where 4^ is the differential flux of the neutrinos given by equation (22), — are the relevant neutrino- nucleon 

charged current differential cross sections, (1 — y) (= E^/ Ei) is the fraction of the neutrino energy transferred to the 
the muon, Vcfi{E^) is the effective volume of the detector, and Tt^"'""^ is the number density of protons in water (= 
Avogadro number). The S-K Collaboration imposed a cut on the upward-going muon path- length of > 7 meters 
in the inner detector which has an effective area of Acff ~ 900 m^ and height « 36.2 m. This 7- meter cut on the 
muon track length can be effectively taken into account by setting Vas = if the effective water-equivalent muon 
range, R^{E^) 5 meters x (E'^/ GeV) , is less than 7 meters, and V^ff = AeS x [i?^(i?^) -I- (36.2—7) meters] 
otherwise [36]. The factor of 1/2 accounts for the fact that only up-going muon events were considered in order to 
avoid the background due to down-going muons produced due to cosmic ray interactions in the Earth's atmosphere. 

The S-K muon events were broadly classified into three categories [49], namely, (i) Fully Contained (FC), (ii) Stop- 
ping (S) and (iii) Through-Going (TG) events. For v^^ energy < 4 GeV the events are predominantly of FC type, 
whereas for energy > 8 GeV the events are predominantly of TG type. Assuming that annihilation of the WIMP of 
mass rriy. produces neutrinos of typical energy ~ (^ — 5)™x' ^® '"^'^ roughly divide the range of our interest into 
three regions according to the resulting muon event types namely, (i) 2 < < 8GcV (FC), (ii) 8 < < 15GcV 
(FC + S) and (iii) 15 GeV <m^{FC + S + TG). 

To set upper limits on the WIMP elastic scattering cross section as a function of WIMP mass for a given annihilation 
channel, we use the following 90% C.L. upper limits [18, 19] on the rates of the upgoing muon events of the three 
different types mentioned above: 



Tp90%C.L. 
'^FC 


~ 13.8 




1 


.^90%C.L. 


~ 1.24 




-1 


Tp90%C.L. 


~ 0.93 


yr^ 


1 



The upper limits on the WIMP-nucleon elastic scattering cross section so derived are then translated into upper limits 
on the branching fractions of various annihilation channels by demanding the consistency of DAMA-compatiblc region 
of the WIMP parameter space with S-K upper limits. These limits are discussed in the next section. 



V. RESULTS AND DISCUSSIONS 



Fig. 1 shows the dependence of the capture rate of WIMPs by the Sun as a function of the WIMP's mass for the 
three halo models specified in Table I. As expected, for a given DM density, the capture rate decreases as WIMP 
mass increases because heavier WIMPs correspond to smaller number density of WIMPs. 

The event rates in the S-K detector as a function of the WIMP mass for the four different WIMP annihilation channels 
are shown in Figure 2 assuming 100% branching ratio for each channel by itself. For each annihilation channel the three 
curves correspond, as indicated, to the three halo models specified in Table I. It is seen that the direct annihilation 
to the lyD channel dominates the event rate, followed by the t'^t~ channel. 

Our main results are contained in Figures 3, 4 and 5, where we show, for the three halo models considered, the 90% 
C.L. upper limits on the WIMP-proton SI and SD elastic cross sections (as a function of WIMP mass) derived from 
the Super-Kamiokande measurements of the up-going muon events from the direction of the Sun [18, 19], for the four 
annihilation channels discussed in the text, assuming 100% branching ratio for each channel by itself. In these Figures, 
we also display, for the respective halo models, the 90% C.L. allowed regions [15] in the WIMP mass vs. WIMP-proton 
elastic cross section plane implied by the DAMA/LIBRA collaboration's claimed annual modulation signal [6], as well 
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FIG. 1: The capture rate as a function of the WIMP mass for the three halo models specified in Table I, and for spin-independent 
(SI: left panel) and spin-dependent (SD: right panel) WIMP-proton interactions. All the curves are for a reference value of the 
WIMP-proton elastic SI or SD cross section of 10~* pb. 
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FIG. 2: The upward going muon event rates in the Super-Kamiokande detector due to neutrinos from WIMP annihilation in the 
Sun as a function of the WIMP mass for the four annihilation channels as indicated, assuming 100% branching ratios for each 
channel by itself, and for spin-independent (SI: left panel) and spin-dependent (SD: right panel) WIMP-proton interactions. 
The three curves for each annihilation channel correspond, as indicated, to the three halo models specified in Table I. All the 
curves are for a reference value of the WIMP-proton elastic SI or SD cross section of 10~* pb. 



as the 90% C.L. upper Hmits [15] on the relevant cross section as a function of the WIMP mass implied by the null 
resuhs from the CRESST-1 [50], CDMS-II-Si [12], CDMS-II-Ge [51] and XENONIO [10] experiments. 

The curves in Figures 3, 4 and 5 allow us to derive upper limits on the branching fractions of the various WIMP 
annihilation channels, from the requirement of consistency of the S-K-implied upper limits on the WIMP-proton 
elastic cross section vifith the DAMA-compatible regions. These upper limits are shovirn in Table II for the three halo 
models discussed in the text. 

Clearly, for the case of spin-independent interaction, there are no constraints on the branching fractions for the bb 
and CC channels since the DAMA-compatible region is already consistent with the S-K upper limit even for 100% 
branching fractions in these channels (the respective curves for the various annihilation channels only move upwards, 
keeping the shape same, as the branching fractions are made smaller). At the same time, for the t~^t~ channel and 
SI interaction, although a 100% branching fraction in this channel allows a part of the DAMA-compatible region to 
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FIG. 3: The 90% C.L. upper limits on the WIMP-proton spin-independent (SI: left panel) and spin-dependent (SD: right 
panel) elastic cross section as a function of WIMP mass derived from the Super-Kamiokande measurements of the up-going 
muon events from the direction of the Sun [18, 19], for the three relevant event types, namely. Fully Contained (FC), Stopping 
(S) and Through Going (TG), as discussed in the text (see Eq. [36]). The thick portions of the curves serve to demarcate the 
approximate ranges where the different event types make dominant contributions to the upper limits. The curves shown 
are for the four annihilation channels, assuming 100% branching ratio for each channel by itself. The 90% C.L. allowed regions 
in the WIMP mass vs. WIMP-proton elastic cross section plane implied by the DAMA/LIBRA experiment's claimed annual 
modulation signal [6] as well as the 90% C.L. upper limits on the cross section as a function of the WIMP mass implied by 
the null results from the CRESST-1 [50], GDMS-ILSi [12], CDMS-II-Ge [51] and XENONIO [10] experiments (solid curves) are 
also shown. All the curves shown are for our halo model Ml (/5dm,q = 0.2GeV/cm=') specified in Table I. 




FIG. 4: Same as Fig. 3, but for the halo model M2 (pdm,© ~ 0.3GeV/cm^) specified in Table I. 



be consistent with the S-K upper limit, consistency of the entire DAMA-compatible region with the S-K upper limit 
requires the branching fraction for this channel to be less than 35-45% depending on the halo model. On the other 
hand, for the channel and SI interaction, there are already strong upper limits (at the level of 25 - 35%) on the 
branching fraction for this channel for consistency of even a part of the DAMA-compatible region with the S-K upper 
limit; and these upper limits become significantly more stringent (by about two orders of magnitude) if the entire 
DAMA-compatible region is required to be consistent with the S-K upper limits. 

The constraints on the branching fractions of various annihilation channels are, however, much more severe in the 
case of spin-dependent interaction; For the quark channels, only parts of the DAMA-compatible region can be made 
consistent with the S-K upper limits, and that only if the branching fractions for these channels are restricted at the 
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(GeV) rrij^ (GeV) 



FIG. 5: Same as Fig. 3, but for the halo model M3 (Pdm,© ~ 0.4 GeV/ cm^) specified in Table I. 





EVENT TYPE 
{m~^ range in GeV) 


UPPER LIMITS ON THE BRANCHING FRACTIONS (in %) (Ml, M2, MS) 


bb 


cc 


T + T- 


vv 


SI 


FC (2.0 - 8.0) 


100, 100, 100 


100, 100, 100 


35,40,45 


0.4, 0.6,0.8 


FC+S (8.0 - 15.0) 


100, 100, 100 


100, 100, 100 


100, 100, 100 


25, 30, 35 


FC+S+TG (15.0 - 20.0) 


100, 100, 100 


100, 100, 100 


100, 100, 100 


25, 30, 35 


SD 


FC (2.0-8.0) 


0.5,0.6,0.7 


0.5,0.6,0.7 


0.05,0.06,0.07 


0.0005,0.0006,0.0007 


FC+S (8.0 - 15.0) 


0.6,0.7,0.8 


0.6,0.7,0.8 


0.14,0.16,0.18 


0.012,0.014,0.016 


FC+S+TG (15.0 - 20.0) 


0.6,0.7,0.8 


0.6,0.7,0.8 


0.14,0.16,0.18 


0.012,0.014,0.016 



TABLE II: Upper limits — derived from Figures 3, 4 and 5 — on the branching fractions for the four annihilation channels, 
from the requirement of consistency of the S-K implied upper limits on the WIMP-proton elastic cross sections with the 
"DAMA-compatible" region of the WIMP mass versus cross section parameter space (within which the annual modulation 
signal observed by the DAM A/LIBRA experiment [6] is compatible with the null results of other DD experiments determined 
within the context of our halo model [15]), for both spin-independent (SI) and spin-dependent (SD) interactions and the three 
halo models specified in Table I. The limits are calculated using the three different upward-going muon event types, namely. 
Fully Contained (FC), Stopping (S) and Through Going (TG). The three consecutive numbers for each annihilation channel 
and muon event type refer to the three different halo models Ml, M2, M3, as indicated. 

level of (0.6 - 0.8)%. On the other hand, for r+r^ and vv channels, parts of the DAMA-compatible regions can be 
consistent vifith S-K upper limits only if their branching fractions are restricted at the level of (0.14 - 0.18)% and 
(0.012 - 0.016)%, respectively, while consistency of the entire DAMA-compatible regions with the S-K upper limits 
requires these fractions to be respectively lower by about a factor of 2.5 (for the r+r^ channel) and a factor of about 
25 (for the vv channel). 

The above small numbers for the upper limits on the branching fractions of the four dominant neutrino producing 
WIMP annihilation channels imply, in the case of spin-dependent WIMP interaction, that the DAMA-allowed region 
of the - cr^p parameter space is essentially ruled out by the S-K upper limit on neutrinos from possible WIMP 
annihilations in the Sun, unless, of course, WIMPs efficiently evaporate from the Sun — which may be the case for 
relatively small mass WIMPs below 4 GeV [36] — or there are other non-standard but dominant WIMP annihilation 
channels that somehow do not eventually produce any significant number of neutrinos while restricting annihilation to 
quark (66, cc) channels to below 0.5% level and r+r^ and vv channels to below 0.05% and 0.0005% level, respectively. 
In the case of spin-independent interaction, however, the DAMA-compatible region of the - ct^^ parameter space 
(or at least a part thereof) remains unaffected by the S-K upper limit if WIMPs annihilate dominantly to quarks 
and/or tau leptons, and annihilation directly to neutrinos is restricted below ^ (25 - 35)% level. At the same time, 
portions of the DAMA-compatible region can be excluded if WIMP annihilation to r+r" occurs at larger than (35 
- 45)% level and/or annihilation to vv occurs at larger than (0.4 - 0.8)% level. These results, based as they are 
on a self-consistent model of the Galaxy's dark matter halo, the parameters of which arc determined by a fit to the 
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rotation curve of the Galaxy, strengthen, at the quahtative level, the earlier conclusion within the SHM that the S-K 
upper limit on the possible flux of neutrinos due to WIMP annihilation in the Sun severely restricts the DAMA region 
of the WIMP mass versus cross section plane, especially in the case of spin-dependent interaction of WIMPs with 
nuclei, although the quantitative restrictions on the WIMP cross section and branching fractions of various WIMP 
annihilation channels obtained here are different (in some cases by more than a factor of 10) from those obtained in 
earlier calculations within the SHM. 



Note added : After the completion of the main calculations of the present work, new results of the S-K collaboration's 
search for upward-going muons ("upmus") due to neutrinos from Sun [20] have appeared. These new results are 
based on a data set consisting of 3109.6 days of data, nearly double the size of the old data set of 1679.6 days used 
in Ref. [18] and in the analysis of this paper so far. Here we consider these new results of Rcf. [20] and the resulting 
changes to the constraints on various WIMP annihilation channels derived above using the earlier S-K results. In 
general, we find that with the new S-K results the upper limits on the branching fractions of various annihilation 
channels become more stringent by a factor of 3 - 4 than those derived above. 



The upmu event categories used in the new S-K paper [20] are somewhat different from those in their earlier work [18]. 
These arc: "stopping" (S), "non-showering through-going" (NSTG), and "showering through-going" (STG); see 
Rcf. [20] for details. For a given WIMP mass, Figure 2 of Ref. [20] allows us to read out the fraction of each 
upmu event type contributing to the total number of events, from the consideration that the typical maximum energy 
of a neutrino produced in the annihilation of a WIMP of mass is ~ 'nixl'^- For low WIMP masses of our interest 
in this paper, rriy^ <, 20GeV (and hence typical neutrino energies ^ lOGeV), the stopping events dominate and con- 
stitute more than 70% of the total number of upmu events, as clear from Figure 2 of Ref. [20]. It is thus expected, 
as indeed we do find from our calculations, that the most stringent upper limits on the branching fractions of various 
WIMP annihilation channels for low WIMP masses come from the observed rate of these Stopping events.^ 



The 90% C.L. Poissonian upper limit on the rate of these Stopping-type upmu events for the new data set of Ref. [20], 
estimated from the total number of this type of upmu events and the number of background upmus due to atmospheric 
neutrinos given in Figure 3 of that reference, is ® ^ 3.27 yr~^. With this, we can calculate, as we did in the analysis 
above, the 90% C.L. upper limits on the WIMP-proton SI and SD clastic cross sections as a function of WIMP mass 
for the new S-K data set of Ref. [20], for the case of 100% branching ratio for each of the four annihilation channels 
considered above. The results, for our best-fit halo model Ml, are shown in Figure 6. 




FIG. 6: Same as Fig. 3, but using the new S-K data from Ref. [20] and considering their "Stopping" upmu events only. 



^ Recall that, for the older data set [18], the most stringent upper limits came from the fully-contained (FC) events; see Table II above. 
* We take the events in the 0-30 degree cone half-angle bin around the Sun to be consistent with the analysis done above for the earlier 
S-K data set. 
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The resulting upper limits on the branching fractions of the four annihilation channels, derived from the requirement 
of consistency of the new S-K-implied upper limits on the WIMP-proton elastic cross sections shown in Figure 6 with 
the "DAMA-compatible" regions, are displayed in Table III. A comparison with the corresponding numbers given in 





Upper limits on the branching fractions (in %) 
from "Stopping" events, with halo model Ml 


bb 


cc 


T + T- 


vv 


SI 


100 


100 


10 


0.11 


SD 


0.12 


0.12 


0.012 


0.00013 



TABLE III: Upper limits on the branching fractions of the four annihilation channels, derived from the requirement of con- 
sistency of the new S-K-implied upper limits on the WIMP-proton elastic cross sections shown in Figure 6 with the "DAMA- 
compatible" regions, for both spin-independent (SI) and spin-dependent (SD) interactions 

Table II shows that the new upper limits on the branching fractions for the relevant annihilation channels are roughly 
a factor of 3-4 more stringent. 

VI. SUMMARY 

Several studies in recent years have brought into focus the possibility that the dark matter may be in the form of a 
relatively light WIMP of mass in the few GeV range. Such light WIMPs with suitably chosen values of the WIMP- 
nucleon SI or SD elastic cross section can be consistent with the annual modulation signal seen in the DAM A/LIBRA 
experiment [6] without conflicting with the null results of other direct-detection experiments. To further probe the 
"DAMA-compatible" regions of the WIMP parameter space — the regions of the WIMP mass versus cross section 
parameter space within which the annual modulation signal observed by the DAM A/LIBRA experiment is compatible 
with the null results of other DD experiments — we have studied in this paper the independent constraints on the 
WIMP-proton SI as well as SD elastic scattering cross section imposed by the upper limit on the neutrino flux 
from WIMP annihilation in the Sun given by the Super-Kamiokande experiment [18, 19]. Assuming approximate 
equilibrium between the capture and annihilation rates of WIMPs in the Sun, we have calculated the 90% C.L. upper 
limits on the WIMP-proton SI and SD elastic cross sections as a function of the WIMP mass for various WIMP 
annihilation channels using the Super-Kamiokande upper limits, and examined the consistency of those limits with 
the 90% C.L. "DAMA-compatible" regions. This we have done within the context of a self-consistent phase-space 
model of the finite-size dark matter halo of the Galaxy, namely, the Truncated Isothermal Model (TIM) [15, 23], in 
which we take into account the mutual gravitational interaction between the dark matter and the observed visible 
matter in a self-consistent manner, with the parameters of the model determined by a fit to the observed rotation 
curve data of the Galaxy. 

We find that the requirement of consistency of the S-K [18, 19] implied upper limits on the WIMP-proton elastic 
cross section as a function of WIMP mass imposes stringent restrictions on the branching fractions of the various 
WIMP annihilation channels. In the case of spin-independent WIMP-proton interaction, the S-K upper limits do 
not place additional restrictions on the DAMA-compatible region of the WIMP parameter space if the WIMPs 
annihilate dominantly to bb and cc, and if direct annihilations to r+r^ and neutrinos are restricted to below ~ 
(35 - 45)% and (0.4 - 0.8)%, respectively. In the case of spin-dependent interactions, on the other hand, the 
restrictions on the branching fractions of various annihilation channels are much more stringent, essentially ruling out 
the DAMA-compatible region of the WIMP parameter space if the relatively low-mass WIMPs under consideration 
annihilate predominantly to any mixture of bb, cc, t^t^, and vD final states. The very latest results from the 
S-K Collaboration [20] put the above conclusions on an even firmer footing by making the above constraints on 
the branching fractions of various WIMP annihilation channels more stringent by roughly a factor of 3-4. Similar 
conclusions were reached earlier [19, 36] within the context of the SHM. The quantitative restrictions on the branching 
fractions for various WIMP annihilation channels obtained here and as given in Table II (and in Table III for the 
latest S-K results [20]) are, however, significantly different from those in the earlier works. 

An important aspect of the Truncated Isothermal model of the Galactic halo used in the present calculation is the 
non-Maxwellian nature of the WIMP velocity distribution in this model, as opposed to the Maxwellian distribution in 
the SHM (see Rcf. [15] for details). This directly affects the WIMP capture rate (and consequently the annihilation 
rate), resulting in significant quantitative differences in the values of the upper limits on the WIMP-proton elastic 
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cross sections (implied by the S-K upper limits on the neutrinos from the Sun) compared to the values in the SHM. 
Similarly, the upper limits on the branching fractions of various possible WIMP annihilation channels (from the 
requirement of compatibility with DAMA results) are also changed. At a qualitative level, however, the general 
conclusion reached earlier [19, 36] within the context of the SHM — that S-K upper limits on neutrinos from the Sun 
severely restrict the DAMA-compatible region of the WIMP parameter space — remains true in the present model 
too, thus adding robustness to this conclusion. 
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